LPTENS-06/22 
LTH-700 
hep-th/0606144 
June 2006 



Chiral family classification of 
fermionic Z<i x Z<i heterotic orbifold models 



Alon E. Faraggi 1 , Costas Kounnas 2 * and John Rizos 3 



1 Dept. of Mathematical Sciences, University of Liverpool, Liverpool L69 1ZL, UK 

2 Lab. Physique Theorique, Ecole Normale Superieure, F-75231 Paris 05, France 
3 Department of Physics, University of Ioannina, GR45110 Ioannina, Greece 



Abstract 

Free fermionic construction of four dimensional string vacua, are related to the 
Z2 x Z2 orbifolds at special points in the moduli space, and yielded the most 
realistic three family string models to date. Using free fermionic construction 
techniques we are able to classify more than 10 10 string vacua by the net 
family and anti-family number. Using a montecarlo technique we find a bell 
shaped distribution that peaks at vanishing net number of chiral families. We 
also observe that ~ 15% of the models have three net chiral families. In 
addition to mirror symmetry we find that the distribution exhibits a symmetry 
under the exchange of (spinor plus anti-spinor) representations with vectorial 
representations . 
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1 Introduction 



The four dimensional superstring vacua based on the free fermionic construction 
[1] are Z 2 x Z 2 toroidal orbifolds at special points in the moduli space [2, 3]. The 
correspondence of the free fermionic point in the moduli space to T-self-dual points 
with maximally enhanced symmetries suggests that symmetry enhancement and self- 
duality play a role in the string vacuum selection. Furthermore, the three generation 
heterotic string models [4] in the free fermionic formulation [1] are the most realistic 
string models constructed to date. The phenomenological appeal, and the theoretical 
considerations, motivate the elaborate study of this class of string compactifications. 

We have therefore embarked in ref. [5] on a complete classification of symmetric 
Z 2 x Z 2 free fermionic orbifold models, according to the chiral content and the four 
dimensional matter gauge group. Thanks to the observation [5] that the twisted mat- 
ter in the models does not depend on the moduli, their chirality classification can be 
carried out at the free fermionic point of the moduli space. Thus, this enables uti- 
lizing free fermionic techniques, which allows an algorithm adaptable to a computer 
program. Resorting to the well known relations in two dimensions between fermionic 
and bosonic currents, one can find the translation of the partition function in the 
bosonic and fermionic representation, and this for any arbitrary point of the moduli 
space. Hence, the free fermionic analysis enables the chirality classification of all 
symmetric, as well as the asymmetric, Z 2 x Z 2 orbifolds. Thus, the free fermionic for- 
malism provides powerful tools for the complete classification of Z 2 x Z 2 perturbative 
string orbifolds. 

The general techniques for carrying out such a classification in the free fermionic 
language were developed in ref. [6] for type II string, and applied in ref. [5] for the 
classification of heterotic chiral Z 2 x Z 2 models. The analysis in ref. [5] was performed 
with respect to a subclass of the models. The Z 2 x Z 2 orbifold of a six dimensional 
compact manifold contains three twisted sectors. In the heterotic string each one 
of these sectors may, or may not, a priori (prior to application of the Generalized 
GSO (GGSO) projections), give rise to spinorial representations. Models that may 
produce, a priori, spinorial representations from all three twisted sectors were dubbed 
S 3 models. This class was classified in ref. [5]. 

It is also possible that the spinorial representations are not present in a given 
twisted plane. Thus, generically we may classify the models in four distinct classes: 
S 3 , S 2 V, SV 2 and V 3 classes of models with spinorial representations arising from 
three, two, one or none of the twisted sectors, respectively. The aim of this work is 
to go beyond the analysis of ref. [5] and complete the chirality classification of the 
Z 2 x Z 2 symmetric orbifolds. 

In this process we find several surprising results. One can envision performing the 
classification by removing or modifying some vectors from the basis set which was 
utilized in ref. [5], in such a way that only two, one or none of the twisted sectors 
produces massless spinorial representations. This method was pursued in ref. [7]. 
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However, it is found that the entire sets of S 3 , S 2 V, SV 2 and V s classes of models 
are produced by working with the basis set of ref. [5] for certain choices of the one- 
loop GGSO projection coefficients (discrete torsions). This result arises from theta 
function identities in the one-loop partition function, which we exhibit in the simplest 
case as an illustration. Hence, these identities allow a complete classification of the 
free fermionic Z 2 x Z 2 orbifold models using a single set of boundary condition basis 
vectors and varying the phases. This enables a systematic analysis of the models 
and the production of algebraic formulas for the main features of the models like 
the number of spinorial, anti-spinorial and vectorial representations. While in the 
past the studies of phenomenologically relevant free fermionic string models has been 
confined to isolated examples, the new methodology allows us to scan a range of 
over 10 16 models, and therefore obtain vital insight into the properties of the entire 
space of Z 2 x Z 2 orbifold vacua. In this paper we present the main outline of the 
analysis and highlights of the results. Further details and results will be reported in 
a forthcoming publication. 



2 Review of the classification method 

In the free fermionic formulation the 4-dimensional heterotic string, in the light- 
cone gauge, is described by 20 left-moving and 44 right-moving two dimensional real 
fermions [1]. A large number of models can be constructed by choosing different 
phases picked up by fermions (/a, A = 1, . . . , 44) when transported along the torus 
non-contractible loops. Each model corresponds to a particular choice of fermion 
phases consistent with modular invariance that can be generated by a set of basis 
vectors Vi, i — 1, . . . , n, 

Vi = {ai(fi),ai(f 2 ),ai(f 3 )) . . .} 
describing the transformation properties of each fermion 

/ A _>_ e iW/A) / A) ,A = 1,...,44. (2.1) 

The basis vectors span a space S which consists of 2 N sectors that give rise to the 
string spectrum. Each sector is given by 

S = '£N i v i , Ni — 0,1 (2.2) 

The spectrum is truncated by a GGSO projection whose action on a string state \S > 
is 



e inVi - Fs \S >=5 S c 



\S>, (2.3) 



where F$ is the fermion number operator and 5s = ±1 is the spacetime spin statistics 



index. Different sets of projection coefficients c 



±1 consistent with modular in- 



variance give rise to different models. Summarizing: a model can be defined uniquely 
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by a set of basis vectors Vi, i 



coefficients c 



, n and a set of 2 N ( N x ^l 2 independent projections 



,i> j. 



The two dimensional free fermions in the light-cone gauge (in the usual notation) 



are: V^;X 



y i u i 1 



l. 



, 6 (real left-moving fermions) and y l ,LU l ,i 



1, . . . , 6 (real 

right-moving fermions), tp A , A = 1, . . . , 5, fj B , B = 1,2,3, (j) a , a = 1, . . . , 8 (complex 
right-moving fermions). The class of models under investigation, is generated by a 
set V of 12 basis vectors 

V = {vi,V2,...,Vi 2 }, 

where 



Vi = 


~- 1 












5 l,..,6^1,..,6^1 A 3^1,..,6 j 






v 2 = 


S 








V2+i = 


ei 


= {y^u^tf,^}, i = 1,...,6, 






Vg = 


h 


= {x 3 \x* 6 ,y 3 \y* 6 \y 3 \f 6 ,f 




.,5 } 


VlO = 


b 2 


= {x 12 ,x 56 ,y 12 ,y 56 \y l \f 6 ,? 




.,6 } 


Vu = 


Zl 








Vl2 = 


Z2 


= {^-.«}. 







(2.4) 



The vectors 1,S generate an ]V = 4 supersymmetric model. The vectors ej,i = 
1, . . . , 6 give rise to all possible symmetric shifts of the six internal fermionized coor- 
dinates (dX 1 = y l u\ dX l = y l u> 1 ), while b\ and b 2 defines the Z 2 x Z 2 orbifold twists. 
The remaining fermions not affected by the action of the previous vectors {S, e\ bi} 
are (ft\i = 1, . . . , 8 which normally give rise to the hidden sector gauge group. The 
vectors Z\ , z 2 divide these eight complex fermions into two sets of four. We stress here 
that the choice of V is the most general set of basis vectors, with symmetric shifts for 
the internal fermions, compatible with a Kac-Moody level one SO (10) embedding. 
Without loss of generality we can fix the associated projection coefficients 



1 




1 




S 




S 




S 




b 2 




S 




= c 




= c 




= c 




= c 






S 


= c 




1 




S 




S 




ei 




b A 






Zn 



= -1, 



leaving 55 independent coefficients, 



,i>3, c 


V 




Z\ 





















ei 




b A 


c 




,c 




,c 






Zn 




b A 




Zn 



i, j = 1, . . .6 , A,B,m,n = 1,2, 



since the remaining projection coefficients are determined by modular invariance [1]. 
Each of the 55 independent coefficients can take two discrete values ±1 and thus a 
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simple counting gives 2 55 (that is approximately 10 16 ' 6 ) distinct models in the class 
of superstring vacua under consideration. 

The vector bosons from the untwisted sector generate an SO (10) x U(l) 3 x 
SO (8) 2 'gauge symmetry. Depending on the choices of the projection coefficients, 
extra gauge bosons may arise from 

x = l + S + Y j e l + z 1 + z 2 = {r? 123 , V> 12345 } 
i=i 

changing the gauge group SO (10) x U(l) — > E$. Additional gauge bosons can arise 
as well from the sectors zi,z 2 and z\ + z 2 and enhance the hidden gauge group 
50(8) 2 -> SO(16) or even 50(8)' -> E 8 . Indeed, as was shown in ref. [5], for 
particular choices of the projection coefficients a variety of gauge groups is obtained. 

The matter spectrum from the untwisted sector is common to all models and 
consists of six vectors of 50(10) and 12 non-Abelian gauge group singlets. The 
chiral spinorial representations arise from the following 48 twisted sectors 

B \i\i\i\ = S + b i + 4e 3 + 4e 4 + l\e 5 + i\e % 

B%&«& = S + b 2 + f iei +e 2 2 e 2 + ele 5 + ele 6 (2.5) 

1 2. a b 

B \t\t\i\ = S + h + %ei + t\e 2 + t\e 3 + t\e A 

where t\ = 0, 1 and b 3 = 1 + S +bi + b 2 + Yli=i £i + J2n=i z n- These states are spinorials 
of 50(10) and one can obtain at maximum one spinorial (16 or 16) per sector and 
thus totally 48 spinorials. Extra non chiral matter, i.e. vectors of 50(10) as well as 
singlets, arise from the B % fi ei ei fi + x , (i — 1,2,3) twisted sectors. 

This construction therefore separates the fixed points of the Z 2 x Z 2 orbifold into 
different sectors. This enables the analysis of the GGSO projection on the spectrum 
from each individual fixed point separately. Hence, depending on the choice of the 
GGSO projection coefficients we can distinguish several possibilities for the spectrum 
from each individual fixed point. For example, in the case of enhancement of the 
50(10) symmetry to Eq each individual fixed point gives rise to spinorial as well as 
vectorial representation of 5O(10) which are embedded in the 27 representation of 
Eq. When E e is broken each fixed point typically will give rise to either spinorial or 
vectorial representation of E 6 . However, there exist also rare situations, depending 
on the choice of GGSO phases, where a fixed point can yield a spinorial as well as 
vectorial representation of 5O(10) without enhancement. The crucial point, however, 
is that the GGSO projections can be written as simple algebraic conditions, and hence 
the classification is amenable to a computerized analysis. 

In ref. [5] we restricted the analysis to the case c ^ = — 1. Prior to GGSO 
projections spinorial representations in this construction can arise from all three 
twisted sectors and we therefore referred to this class as 5 3 models. This is somewhat 
of a misnomer as we discuss below. To produce models with spinorial representations 
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arising only from one or two of the twisted one can contemplate modifying the basis 
vectors. For example, removing z 2 from the set will entail that the third twisted place 
produces only massive states. Hence, this would correspond to the models dubbed 
as S 2 V. Similarly, modifying the Z 2 x Z 2 basis vectors b\ and b 2 in such a way that 
they produce vectorial rather than spinorial representation and removing both z\ and 
z 2 from the base would entail that only vectorial representations are generated. An 
analysis along this lines was followed in ref . [7] . However, as a result of Jacobi theta 
function identities it turns out that the analysis can be carried entirely with the basis 
(2.4) and just modifying the phases. To illustrate this correspondence we consider 
the simplest possibility given by the set {1, S, x}. The partition function is: 



{91 - e$ - etmef + e%ef + m*m + <% + m 



(2.6) 



The gauge group is SO (28) x E 8 . Now consider the set {1, S, x, zi}. There are now two 
consistent choices of the coefficient c x — ±1. The choice c x = +1 produces the 

.21 J L z lJ 

SO (28) x E 8 gauge group, while the choice c * = -1 produces an SO(20) x 50(24) 
gauge group. Indeed, the one-loop partition function as function of c 

' Jtf 4 #4 M\ fa6nW\ f 7)4 7)8 , 7)47)8 , 7)47)8, 
7, \ y 3 _t/ 2 _y 4/ \ y 3 y 3 / \ y 3 y 2 + ^3^3 + y 3 y 4 + 



becomes 



Mp8 _j_ Mn8 _j_ Mn8, 

u 4 u 2 + y A y s + u A v A + 



Q\Q\ + ^2^3 C 



/ 4 i/ 3 
^1 



7 4 u 4 



m } + 



Plus two additional groups of terms with the permutation of 3 , 6 2 and 6±. Fixing 
c ^ = +1 and using the Jacobi identity 9\ — 9 2 — 9\ = reproduces the partition 
function of the set {1, S, x}. This simple example illustrates how we may obtain vacua 
from an enlarged basis set {1, S, x, Z\\ identical to those obtained from a reduced 
basis set {1, S, x} for an appropriate choice of the GGSO coefficient c 
primary feature which is exploited in our classification. 



and is the 



3 Counting the twisted matter spectrum 

The counting of spinorials can proceed as follows. For each 50(10) spinorial B l pqrs 
in (2.5) we write down the associated projector Pp qrs = 0, 1. The detailed expressions 
for the 48 projectors are 



P (i) 




1 - c 



1 -c 



B 



B, 



e 2 
(i) 

z 2 

(1) 



6 



3(2) 



3 (3) 



1 

4 
1 
4 
1 
4 
1 
4 



— c 



— c 



e 3 

o(2) 
o(2) 

e 5 

o(3) 

.-^l-MsAl. 

Zl 

o(3) 



— c 



— c 



e 4 

o(2) 
\Pl\i7.lht%. 

Z2 
o(2) 

e 6 



(3-1) 



5 



(3) 

^2 
o(3) 



For the surviving spinorial {Pl qr s = 1) the chirality (16 or 16) is determined from 



the associated chirality coefficient X, 



pqrs 



±1, where 



X, 



(i) 



(2) 



= — c 



= — c 



h + (1 - 4)e 5 + (1 

o(l) 

h + (1 - 4)e 5 + (1 



X (3) - -r 



6)e 6 ' 

\ + (1 - £ 3 )e 3 + (1 - 4)e 4 



0(2) 



(3.2) 



o(3) 



Using the above results, we can easily calculate the number of spinorials/antispinorial 
per sector 



1 ± 

E P9 rs p(i) j - 1 9 Q 



(3.3) 



The counting of SOiVS) vectorials can proceed in a similar way. For each vectorial 
generating sector B pqrs + x the associated projector P pqrs is obtained from (3.1) using 
the replacement B pqrs — > B pqrs + x. Since there is no chirality in this case the number 
of vectorials per sector is just the sum of the projectors 



y{i) _ V P(* 



) 

pqrs 



(3.4) 



pqrs 



The total vectorial (F) and net spinorial (S) numbers are 



v = J2 



(3.5) 



i=l 



and 



(0 

+ 



s 



00 



(3.6) 
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The mixed projection coefficients entering the above formulas can be decomposed 



in terms of the independent phases c 



> j. After some algebra we come to the 



conclusion that for the counting of the spinorial/antispinorial and vectorial SO (10) 



states the phases c 



= 1,...,6, 



,A = 1,...,2, 



as c 
tota 



are not relevant. Moreover the phase 



1, . . . , 2 as well 
is related to the 

chirality flip. This leaves a set of 40 independent phases which is still too large 
for a manageable computer analysis. We therefore resort to a Monte-Carlo analysis 
that generates random choices of phases. The complete classification is currently 
underway and will be reported in a future publication. In this sense our results are 
based on a statistical polling, with a sample that correspond to some 10 10 vacua. We 
have checked however that increasing the size of the sample by 10% does not alter 
our results. Therefore our results, while in some sense empirical, are expected to hold 
for the entire space of vacua. 



4 The observable and hidden gauge groups 

An important step in the analysis is the determination of the four dimensional 
gauge group. In particular, it is important to determine the component of the 
gauge group which is identified with the observable gauge group. Observable mat- 
ter then consist of the states that are charged with respect to this gauge group. 
In our construction the sectors that produce space-time vector bosons include: 
G = {0, zi, Z2, z\ + Z2,x}. The sector gauge bosons produce the gauge group 
SO (10) x U(l) 3 x SO(8) 2 . Depending on the choice of the GGSO projection coeffi- 
cients the four dimensional gauge group is enhanced. Since the sectors that produce 
space-time vector bosons do not break any supersymmetries the classification of the 
four dimensional gauge group can be done at the level of the N = A vacuum. The ad- 
dition of the supersymmetry breaking sectors b\ and 6 2 then breaks the N = 4 gauge 
group to the unbroken gauge group at the N = 1 level. Some of the possibilities at 
the N — 1 level are listed in table 1. 

Several features are noted from the list. The first five groups are descendants of 
SO (32), whereas the later five are descendants of E$x E$. This gross division is con- 
trolled by the phase c ^ . An important point to note, relevant for our classification, 
is the occurrence of models with several 5*0(10) group factors. This arises because 
of the enhancement of the 0-sector 5*0(8) x U(l) group factors to 5O(10). In our 
analysis we define the observable gauge group to be 5*0(10) and the chiral matter 
should be charged under that group. The question arises as to which 5O(10) should 
be identified as the observable gauge group and the subsequent classification of the 
chiral and vectorial matter states with respect to that group. In our analysis here 
we identity the observable 5O(10) symmetry as the one which is generated by the 
world-sheet fermions ip 12345 . Investigations of other possibilities are left for future 
work, although our experience with the construction of quasi-realistic string models 
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Gauge group 
£0(10) x £0(18) x U(lf 
£0(10) x £0(9)" x U(lf 
_ £O(10) 2 x £0(9) x U(lf 
£O(10) 3 x £7(1) 
£0(26) x U(lf 
E 6 x U(lf x E 8 
E e x U(lf x £0(16) 
E 6 x U(lf x £0(8) x £0(8) 
£0(10) x U(lf x E 8 
£0(10) x U(lf x £0(16") 

Table 1: Typical enhanced gauge groups for a generic model generated by the basis 
(24)- 

suggests that the spectrum with respect to other £0(10) group factors is vectorial 
rather than spinorial. 

We further comment that in our classification the observable GUT gauge group is 
always £0(10). The conditions for enhancement of the gauge group are given in ref. 
[5], eqs. (4.3-4.11). These conditions are incorporated into our classification routine 
and cases in which the £0(10) symmetry is enhanced are rejected. 

We also note that our classification is with respect to the chiral content of the 
models, which in the free fermionic models arises from the twisted sectors. The break- 
ing of the GUT £0(10) symmetry to a subgroup, which can be further broken to the 
Standard Model gauge group, is achieved in free fermionic models with an additional 
boundary condition basis. Free fermionic models with quasi-realistic gauge group and 
chiral family content were produced in the past [4]. The interest in this paper is in 
the global properties of a large class of compactifications to which the quasi-realistic 
free fermionic models belong, but not in producing quasi-realistic spectrum. The 
classification of free fermionic models with broken £0(10) GUT symmetry will be 
pursued in future work. 

5 Results 

The results of the random search are exhibited in figures 5 and 2. The first figure 
shows the percentage of models with a net number of chiral families. The second 
figure exhibits the total number of vacua on a logarithmic scale over the net number 
of chiral families. In our sample the peak of the distribution is for a vanishing net 
number of chiral families. About 15% of the models contain a net number of three 
chiral or anti-chiral families. By increasing the sample size by 10% we note that 
these results are not modified. 
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Figure 1: Distribution of a Monte Carlo generation of models sampling some 10 10 
vacua. 

In figure 2 we present scatter plot of the logarithm of the number of models 
versus the net number of chiral families. The vertical spread arises from variations 
in the number of vectorial representations in the models. The peak is for vanishing 
net number of chiral families, and we note that the distribution is symmetric about 
that point. The plot has a bell shape which recedes for increasing number of chiral 
families. For net family number above 24 the number of occurrences is small and 
are not shown on the plot. This reflects that increasing number of chiral families 
requires that the configurations of the GGSO phases attain higher symmetry and 
consequently the number of possibilities decreases. Curiously we note that models 
with some net number of chiral families do not appear on the plot. For example 
models with 7, 9, 11, 13, 14, 15, 17, 18, 19 net number of chiral families do not appear 
in our sample. Thus, these are either very rare or forbidden altogether in symmetric 
Z 2 x Z 2 orbifolds. Additional plots and analysis of the data will be presented in a 
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Figure 2: Scatter plot of the logarithm of the number of models versus the net number 
of chiral families. 

forthcoming publication. 

The symmetry about the vanishing number of chiral families is in accordance 
with mirror symmetry. We note, however, an additional symmetry in the distribu- 
tion under exchange of vectorial, and spinorial plus anti-spinorial, representations. 
The symmetry states that for a model with a given total number of spinorial plus 
anti-spinorial representations there exist a corresponding model in which the spino- 
rial plus anti-spinorial representations are exchanged with vectorial representations. 
Ultimately, in the free fermionic language this symmetry reflects a symmetry under 
a discrete exchange of some GGSO projection coefficients. However, it may have 
interesting implications in terms of the underlying geometrical data. 

In conclusion, the quasi-realistic free fermionic models are among the most re- 
alistic string models constructed to date. While past exploration of these models 
consisted of the study of single examples, in this paper we embarked on the inves- 
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tigation of the properties of the whole space of vacua in this class. Future studies 
will incorporate into the analysis further properties of the realistic models, like the 
assignment of asymmetric boundary conditions. While the task is still horrendous in 
terms of the shear number of vacua, the existence of string models in this class that 
come close to describing reality, and precisely where one would most expect to find 
them, gives ample reason to suggest that we are on the right track. 
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